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We suggest a method for creating bright matter solitons by loading a BEC of atoms in a driven 
tilted optical lattice. It is shown that one can realize the self- focussing regime for the wave-packet 
dynamics by properly adjusting the phase of the driving field with respect to the phase of Bloch 
oscillations. If atom-atom interactions are larger than some critical value g m in, this self- focussing 
regime is followed by the formation of bright solitons. Increasing the interactions above another 
critical value g m ax makes this process unstable. Instead of soliton formation one now meets the 
phenomenon of incoherent self-trapping. In this regime a fraction of atoms is trapped in incoherent 
localized wave-packets, while the remaining atoms spread ballistically. 



Introduction. Non-spreading localized wave-packets, 
solitons, are a paradigm of nonlinear wave dynamics and 
are encountered in many different fields, such as physics, 
biology, oceanography, and telecommunication [1]. Re- 
cently much attention has been paid to solitons with a 
Bose-Einstein condensate (BEC) of cold atoms [3-01 , de- 
scribed in the mean-field approximation by the nonlinear 
Schrodinder equation. To generate solitons this equation 
should be self-focussing, which requires either attractive 
inter- atomic interactions or negative mass. The latter 
case can be achieved by loading a BEC into an optical 
lattice, where the effective atomic mass is negative near 
the edge of the Brillouin zone. This idea was realized in 
experiment 0] where a wave-packet of cold atoms was 
moved to the edge of the zone by accelerating the op- 
tical lattice for a given time. Adjusting the number of 
atoms in the wave-packet and its width the authors of the 
cited experiment were able to create a stationary soliton, 
which showed no sign of dispersion for at least 60ms. 

In this work we suggest a different method for creating 
bright solitons by loading a BEC of atoms in a driven 
tilted optical lattice. This method relies on the well- 
known result that the quasi-energy spectrum of the sys- 
tem has band structure if the driving frequency coincides 
with the Bloch frequency (see for example review |8[ and 
references therein). Changing the phase of the driving 
field with respect to the phase of Bloch oscillations it is 
possible to realize both defocussing and focussing regimes 
of the nonlinear Schrodinder equation. We demonstrate 
that the suggested method is robust and does not re- 
quire the simultaneous adjusting of wave-packet width 
and atom number. Moreover, it allows us to generate 
moving solitons, which is a challenging problem from the 
view point of laboratory experiments. 

The model Considering deep enough optical lattices, the 
BEC dynamics in a tilted driven lattice obeys the follow- 
ing equation 

ifiai = -^(az + i+az_i)+d[F-F w cos(cjt+0)]/ai+p|ai| 2 a^ , 

(1) 

where a\ is the BEC complex amplitude in /th potential 
well, J the hopping matrix element, d the lattice period, 



F and F u magnitudes of a static and AC fields, uj the 
driving frequency, and <j) an arbitrary phase. For van- 
ishing atom-atom interactions system (pQ) has been dis- 
cussed in numerous papers, including the experimental 
works [9I4TH. For g ^ and arbitrary values of the 
driving frequency system ffl has been studied in the 
recent experimental work [12| and the theoretical work 
[l3|. Here we focus on the particular case where the 
driving frequency uj coincides with the Bloch frequency 
uob = dF/h. As initial conditions we consider a coher- 
ent wave-packet of a finite width and, to be specific, we 
assume ai(t = 0) ~ exp(— I 2 /2o~q). 

One gets a useful insight into the BEC dynamics by 
changing from the Wannier basis to the basis of Wannier- 
Stark states, which corresponds to the substitution 

c m = J2^m-i(J/dF)a l . (2) 
1 

Assuming additionally F u <C F and using the rotating- 
wave approximation, the equation with new amplitudes 
c m takes the form 

ihdm = -^(c m+:L e^ + c m _ie"^) 
^2 Im,m 1 ,m 2 ,m 3 c mi c m2 c rn3 5(m+m 1 -m2-m 3 ) , 

mi ,'Ti2 ,7^3 

(3) 

where J = JF 0J /2F and the kernel Im,m 1 ,m 2 ,m 3 — 

Y,n^n-rn(z)Jn-m 1 (z)Jn-rn 2 (z)Jn-rri3( z )- ( Here and 

above J n (z) are Bessel functions of the first kind and 
z = J/dF.) Since we have assumed uj = ujb, the param- 
eter (j> in Eq. (|3j) takes into account the phase difference 
between the Bloch and field oscillations. 

Dynamics for dF ^> J. We begin our analysis of 
Eq. (j3j) by considering the case of a strong static field, 
dF ^> J. In this limiting case the Wannier- Stark states 
coincide with the Wannier states and Eq. (j3j) reduces to 
the celebrated Discrete Nonlinear Schrodinger Equation 

ihc m = -^(c m+1 e i(f) + c m _ie" i0 ) + #|c m | 2 c m . (4) 
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This equation can be considered as a discretization of the 
continuous nonlinear Schrodinger equation 
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Let us discuss the system ((4]) in terms of its continuous 
counterpart (j5j). The dynamical regimes of (j5j) depend 
on the value of the parameter <p. The most interesting 
cases are </> = (defocussing of the initial wave-packet), 
(j) = ±7r/2 (translation), and <p> = ir (self- focussing). It 
is easy to check numerically that under the condition 
do ^> 1 these regimes are present in the discrete equation 
dl]) as well. From now on we will focus on the case <fi = ir 
where the effective hopping matrix element is negative 
and Eq. (|5]) is self-focussing. As known in the ID nonlin- 
ear Schrodinger equation the transient focussing regime 
is followed by formation of several solitons moving with 
different velocities [14[. In the co- moving frame each of 
isolated solitons is approximated by the one-parameter 
function, 



f(x;X) 
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(6) 



which solves the stationary nonlinear Schrodinger equa- 
tion. The parameter A in (|6j), which defines the soliton 
width, is proportional to the number of atoms captured 
by the soliton and the nonlinearity g. Namely, assuming 
the normalization of the wave function ijj(x,t) to unity, 
we have A = ag/2Jd, where a is the relative number 
of atoms in a given soliton. Simulating the dynamics 
of a single (moving) soliton one finds the function (j6j) 
with discrete x = md to be an approximate solution of 
dl]) until A « 0.1, where the soliton width is about 20 
lattice periods. For a larger A the discreetness of x be- 
comes important and one speaks about a discrete soliton 
or breather [HI, [l6|. The discrete soliton may involve 
only 3-4 sites and is pinned in the lattice. 

We simulated the DNLSE (j4]) for a wide initial wave 
packet, searching for stable soliton-like structures, and 
indeed, observed them in a certain interval of nonlinear- 
ity gmin < 9 < Qmax- An example is given in Fig. [fl 
where g = g/J = 2. (Since J in Eq. (j4]) can be set to 
unity by rescaling the time, only the ratio g/J matters.) 
In Fig. [T]one sees two solitons of the form ([6]) moving in 
opposite directions with constant velocities, and one cen- 
tral soliton, pinned at the lattice origin. This central soli- 
ton, which captures the largest part of atoms, is actually 
on the border of validity of the continuous approximation 
and carries some features of discrete solitons. Increasing 
the interaction constant g more atoms become captured 
in the central soliton-like structure and it becomes a dis- 
crete breather with approximately three sites involved in 
the dynamics (see upper panel in Fig. [2]). The remaining 
atoms are emitted in the form of plane waves and tiny 
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FIG. 1: Lower panel: Dynamics of J4j) for <p = 7r, gq = 20, 
J = 1, and g — 2. The time is measured in units of the 
tunneling period T — 2nTi/ J, which defines the characteristic 
time scale of the system dynamics. The upper limit of the 
color axis for the gray-scaled encoding of site populations is 
set to 0.04. Upper panels: Site populations P m = |c m | 2 at the 
end of numerical simulations for N = oo (left) and N = 10 6 
(right). 



solitons. However, these solitons appear to be sensitive 
to the interaction-induced decoherence process present in 
the system and, hence, are of little interest. 

Next we discuss the destructive effect of the 
interaction-induced decoherence in more detail. To 
mimic the quantum decoherence within the framework 
of the DNLSE the solution of (pQ) should be averaged 
over an ensemble of initial conditions. This is uniquely 
defined by the initial quantum many-body state of the 
system [13]. In the case of the BEC initial state being 
considered at present, this ensemble can be approximated 
by the following ensemble, 



cJ = e*VM 2 + ej, ef = \ Cl \ 2 /N, 6? = l/W\ Cl \ 2 , 

(7) 

where e\ and Q\ are normally distributed random variables 
and N is the total number of atoms. Clearly, the ensem- 
ble ([7]) takes into account quantum fluctuations of atom 
number in any given well of the optical lattice. For large 
N and providing the solution of the DNLSE is stable with 
respect to a small variation in initial conditions, this ad- 
ditional averaging procedure may be omitted. However, 
if the solution is not stable this averaging procedure is ab- 
solutely necessary. Firstly, because it removes artifacts of 
the mean- field discription [13, [H| and secondly, because 
it makes transparent the decoherence of the BEC in the 
course of time [19|. 

The upper panels in Fig. [1] compare the single-run so- 
lution of DNLSE with that averaged over the ensemble 
0, where we set TV = 10 6 . It is seen that the averaging 
procedure only slightly increases the soliton widths; in all 
other aspects the time evolution of site populations fol- 
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FIG. 2: The same as in Fig. \T\ yet g = 4. The upper and 
lower panels refer to N = oo and N = 10 6 , respectively. The 
central breather in the upper panel captures approx. 40% of 
the atoms, the largest side solitons approx. 5% each. 



low that depicted in the lower panel in Fig. [TJ However, 
this result holds only if g < g m ax, where the system dy- 
namics are stable (gmax ~ 2 for the chosen parameters). 
A comparison of the upper and lower panels in Fig. [2j 
where g > g m ax and the system dynamics are unstable, 
indicates that the decoherence process destroys the small 
side solitons. Simultaneously it transforms the large cen- 
tral breather into an incoherent packet. For this reason 
we avoid using the term 'discrete soliton' or 'breather', 
which is reserved in physical literature for a coherent lo- 
calized excitation. Instead we shall use the term 'self- 
trapping', which may be used for both coherent and inco- 
herent localized wave-packets. It is also worth mention- 
ing that the discussed (incoherent) self-trapping princi- 
pally differs from a (coherent) self-trapping in the defo- 
cussing regime (positive hopping matrix element) studied 
in Refs. [20|, Hl|. We come back to the problem of inco- 
herent self-trapping in the next section. 

Dynamics for dF < J. We also studied the forma- 
tion of bright solitons in the driven tilted lattice (pQ) in 
the parameter region, where it cannot be reduced to the 
standard DNLSE. In particular, when relaxing the con- 
dition dF ^> J the sum in the nonlinear term in Eq. (|3]) 
contains many terms; thus, one cannot appeal to the non- 
linear Schrodinger equation. Nevertheless, we find the 
dynamics of system ([3]) to be similar to that of system 
Namely, for <fi = one observes defocusing of the ini- 
tial wave packet, translation dynamics for <p = 7r/2, and 
self- focusing for <j) = tt. Moreover, the focussing regime 
is followed by formation of a number of solitons with 
the characteristic shape given in Eq. (|6j). An example is 
given in the lower panel in Fig. [3j where J/dF = 4 and 
the sum in Eq. ([2]) effectively contains 10 terms. Note 
that this figure depicts populations of the Wannier- Stark 
states. Populations of the lattice sites show additional 



FIG. 3: Dynamics of the original system (pQ) for <j> = 7r, a = 20, 
N 10 6 , J = 2, dF = 0.5, Fu 0.2F, and g 0.8 in 
the Wannier (upper panel) and Wannier-Stark (lower panel) 
basis. The chosen parameters correspond to J = 0.2 and 
g = 4. The time is measured in units of the effective tunneling 
period T = 2nh/J 



oscillations with the Bloch frequency (see upper panel). 

The result depicted in Fig. [3] refers to moderate non- 
linearity g = 0.8 or g = g/ J = 4, where the system dy- 
namics are stable. Larger atom-atom interactions lead to 
unstable dynamics and cause a rapid decoherence of the 
initial BEC state. Then, instead of soliton formation, one 
meets the phenomenon of temporal self-trapping. This 
phenomenon is exemplified in the lower panel in Fig. 3J 
which shows evolution of the site populations for g = 100. 
It is seen that atoms are temporally trapped by the non- 
linearity. A convenient characteristic of this process is 
the survival probability S(t), which we define as the rel- 
ative number of atoms staying in the region of support 
of the initial wave-packet (\l\ < 50 for <jo = 20). The 
function S(t) is depicted by the solid line in the upper 
panel in Fig.HJ It should be noted that S(t) goes to zero 
for t — » oo. This constitutes the main difference between 
the self-trapping in the DNLSE regime (dF ^> J) and in 
the regime where dF < J; in the former case S(t) ap- 
proaches a constant value defined by the relative number 
of permanently trapped atoms. 

Finally, we would like to stress one more time that the 
trapped atoms form an incoherent wave-packet, where 
the single-particle density matrix is close to a diagonal 
matrix. In fact, for a systematic study of the incoherent 
self-trapping it is desirable to have an incoherent wave- 
packet from the very beginning. For the sake of com- 
parison the dash-dotted and dashed lines in the upper 
panel in Fig. 2] depict S(t) for completely incoherent ini- 
tial conditions and g = and g = 100, respectively (22j . 
In the former case the self-trapping is absent and the 
atoms spread ballistically [13]. In the latter case, decay 
of the survival probability practically coincides with that 
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FIG. 4: Lower panel: Dynamics of the original system (pQ) in 
the Wannier basis for = tt, a = 20, N = 10 6 , J = 2, F = 
0.5, Fcj = 0.1F, and g = 10. These parameters correspond to 
J = 0.1 and g = 100. Upper panel: Survival probability for 
the coherent initial wave-packet (0) and g = 100 (solid line), 
for the incoherent initial wave-packet and g — 100 (dashed 
line), and for the incoherent initial wave-packet and g — 
(dash-dotted line). 

depicted by the solid line if we take into account the time 
needed to decohere the initial BEC state. 

Conclusions. We have analyzed the dynamics of a BEC 
of cold atoms with repulsive interactions in driven tilted 
optical lattices. It is shown that one can easily realize the 
focussing regime for the wave-packet evolution by prop- 
erly adjusting the phase of the driving field with respect 
to the phase of Bloch oscillations. If the macroscopic 
interaction constant g is larger than some critical value 
gmin, this focussing regime is followed by the formation 



of several bright solitons. We also checked that these 
solitons are stable against the interaction-induced deco- 
herence process present in the system. Further increase 
of the interaction constant above another critical value 
Qmax greatly enhances the decoherence process and, in- 
stead of soliton formation, one meets the phenomenon 
of incoherent self-trapping. In this regime a fraction of 
atoms is permanently or temporally trapped in an inco- 
herent localized wave-packet, while the remaining atoms 
spread ballistically with the velocities defined by the ef- 
fective hopping matrix element J. 

The discussed incoherent self-trapping has been re- 
cently observed in the laboratory experiment [12| . In the 
cited experiment the BEC of N = 1.2 x 10 5 cesium atoms 
in vertical quasi one-dimentional lattice of the depth 3 
recoil energies was subject to the sum of gravitational 
field and levitation force due to a magnetic field gradi- 
ent, resulting in the Bloch frequency ujb — 27r x 98 Hz. 
Modulating the magnetic field with the same frequency, 
the authors of Ref. [12| realized the model (pQ), where the 
ID macroscopic interaction constant g (which is propor- 
tional to the s- wave scattering length a s ) was varied in a 
wide interval by means of the Feshbach resonance. The 
regime of incoherent self-trapping was clearly observed at 
a s = 336ao (ao is Bohr radii), where the rate of spread- 
ing of the atomic cloud was essentially suppressed in con- 
strast to the case a s = llao- Simulating the dynamics of 
system (pQ) for the specified parameters we concluded that 
to enter the regime of stable bright solitons, one should 
increase the number of atoms to N ~ 10 6 and further 
decrease the ID interaction constant g. 
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